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ABSTRACT 

2 

We  consider  the  one  dimensional  porous  media  problem  *  (u  )xx, 
u(0,x)  -  Uq ( x ) ,  with  initial  data  uQ  compactly  supported.  We  deduce  a 
coordinate  transformation  x  ”  X(t.p)  which  renders  stationary  the  free 
boundary  separating  {(t,x):  u(t,x)  >  0}  and  {(t,x):  u(t,x)  >0}.  We 
describe  a  simple  difference  scheme  based  on  this  transformation  and  show 
that,  even  with  a  fairly  crude  mesh,  the  free  boundary  is  tracked  quite 
accurately. 
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SIGNIFICANCE  AND  EXPLANATION 

[1  ^ 

We  consider  the  one-dimensional  porous  media  problem  -  (u»xx, 
u(0,x)  ■  Uq(x)  for  initial  data  ^  which  vanishes  outside  a  finite 
interval.  Solutions  of  this  problem  contain  a  free  boundary  separating 
points  (t,x)  with  u(t,x)  >  0  from  those  with  u(t,x)  ”  0.  We  deduce  a 
coordinate  transformation  x  -  X(t,p)  which  renders  this  free  boundary 
stationary  in  time.  We  describe  a  simple  difference  scheme  based  on  this 
transformation  and  show  that,  even  with  a  fairly  crude  mesh,  the  free  boundary 
is  tracked  quite  accurately. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
susmMry  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


A  COORDINATE  TRANSFORMATION  FOR  THE  POROUS 
MEDIA  EQUATION  THAT  RENDERS  THE  FREE-BOUNDARY  STATIONARY 


Morton  E.  Gurtin,  Richard  C.  MacCany ,  and  Eduardo  A.  Socolovsky 
1.  Introduction . 

In  this  paper  we  study  degenerate  diffusion  problems  in  which  free 
boundaries  occur ,  our  major  objective  being  the  development  of  numerical 
procedures  which  effectively  track  these  boundaries*  For  convenience,  we 
introduce  our  ideas  in  terms  of  the  one-dimensional  porous  media  problem; 
generalizations  are  given  at  the  end  of  Section  2  and  in  Section  4. 

Consider  the  problem  of  determining  u(t,x),  t  >  0,  x  e  R,  such  that 


ut  -  (u2) 
u(0,x)  « 


XX 


u0(x) 


$ 


(P) 


with  initial  data  uQ  supported  on  a  finite  interval: 

uQ(x)  >0,  -a  <  x  <  a  , 


(1.1) 


uQ(x)  »  0  otherwise 

This  is  (a  special  case  of)  the  porous  media  problem  (cf«,  e.g.,  [1]).  As  is 
known  ( [2  -  6] ) :  there  exists  a  unique  weak  solution  u;  at  each  t  the 
support  of  u  is  a  finite  interval 


C_(t)  <  x  <  c+(t)  » 

inside  its  support  u  is  smooth,  but  across  C^(t),  ux  is  generally 
discontinuous; 

5+(t)  -  -2ux(t,C±(t))  .  (1.2) 
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There  exist  numerical  procedures  for  Problem  (P)  (cf.,  e.g.,  [7  -  9]), 
all  of  which  seem  to  encounter  some  difficulty  in  tracking  the  free  boundaries 
x  -  C^(t).  The  main  feature  of  our  method  is  the  determination  of  a  family  of 
curves 


x  =  X(t,p) ,  X(0,p)  -  p 

along  which  the  free  boundary  propagates: 

?±<t)  -  X(t,+a)  . 

Writing  (  P)  1  in  the  form  of  a  mass  balance  law 

ut  +  (uv)x  -  0 

with  u(t,x)  a  "density"  and 


(1.3) 


(1.4) 


v  -  -2ux  (1.5) 

a  "velocity"  (cf.,  e.g.,  [10]),  we  see,  using  (1.2),  that  the  free  boundary 
propagates  with  the  velocity  v(t,x)  of  the  medium;  hence  property  (1.3)  will 
follow  provided  we  take  X(t,p)  to  be  the  solution  of 

Xt(t,p)  -  v(t,X(t,p) ) ,  X(0,p)  -  p  .  (1.6) 

Within  this  context  x  *  X(t,p)  represents  the  motion  of  the  mediua  with 
material  points  labeled  by  their  positions  p  at  t  *  0. 

Using  these  ideas,  we  are  able  to  reduce  ( P)  to  the  following  initial- 
value  problem  for  X: 


XX-  2u.X  -  2u' X  , 
p  t  0  pp  Op 

X(0,p)  -  p  . 


(P*) 


ID  m 

(Here  and  in  what  follows  X  -  (X  )  .)  The  free  boundary  is  then  given  by 

P  P 


U  -  X_1uft  , 
P  0 


(1.3),  while  U(t,p)  -  u(t,X(t,p))  satisfies 

(1.7) 

a  relation  which  expresses  balance  of  mass  in  material  (Lagrangian) 
coordinates.  Our  procedure  for  solving  (  )  consists  in  solving  (P  )  on  the 
fixed  interval  -a  <  p  <  a  of  support  of  Uq. 
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In  Section  2  we  establish  a  uniqueness  theorem  for  ( P*) ,  and  we  show  that 

given  a  sufficiently  regular  solution  X  of  this  problem,  (1.7)  generates ,  at 

least  locally  in  tine,  a  weak  solution  u  of  our  original  problem  (P). 

The  determination  of  u( t,x)  from  0(t,p)  requires  that  X(t,p) ,  as  a 

function  of  p,  be  invertible  at  each  t ,  a  condition  related  to  the  non- 

vanishing  of  X .  We  show,  in  Section  2,  that  X  >  1  for  all  time  whenever 
P  p 

u”  <  0  on  its  support,  and  that  Xp{t,a),  say,  tends  to  zero  in  a  finite 

time  T  whenever  u^(a)  ”  0,  ujj(a)  >  0.  We  show  further  that  under  the 

latter  two  conditions  X(t,a)  =  a  for  0  <  t  <  T,  so  that  the  free  boundary 

is  vertical  until  t  »  T  (cf.  [6,  11)).  We  also  establish  a  growth  estimate 
2 

for  the  L  (-a, a)  norm  of  X  (t,*). 

P 

* 

In  Section  3  we  describe  a  simple  difference  scheme  for  Problem  (P  )  and 
give  some  calculations  which  demonstrate  the  utility  of  our  procedure;  in 
particular,  we  show  that  even  with  a  fairly  crude  mesh  the  free  boundary  is 
tracked  quite  accurately. 

While  our  paper  is  devoted  to  the  one- dimensional  porous  media  problem, 
our  method  seems  to  have  considerable  generality;  in  Section  4  we  derive  the 
analog  of  (P  )  for  the  porous  media  problem  in  sf1 ;  in  a  future  paper  we  will 
discuss  applications  to  more  general  equations  and  to  Stefan  problems. 


We  first  proceed  formally.  Let  u  be  a  solution  of  (P)  with  initial 
data  Ug  subject  to  (1.1),  let  X(t,p)  be  the  solution  of  the  initi al-value 


problem  (1.6),  and  define 


D(t,p>  ■  u(t,X(t,p)) 

or  more  generally,  for  any  function  f(t,x), 

f#(t,p)  -  f(t,X(t,p)) 
Then  by  (1.4),  (1.6),  and  (2.1), 

Ut  "  <ut  +  uxv)#  "  -(uvx)#»  V 

so  that 


<v\ 


(2.1) 


(2.2) 


(0«p)t  -  0t2p  ♦  t«pt  -  0 

But  by  (1.6)2, 


(2.3) 


n(o,p)  -  u(0,p)  ■  u0(p),  xp(o,p)  -1  i 

hence 

0  -  X~'un  .  (2.4) 

p  0 

Further,  by  (2.1), 

°p  ■  <v%  • 

and  (1.5),  (1.6)  yield 

X  -  -20  x’1  .  (2.5) 

t  P  P 

Equations  (2.4)  and  (2.5)  form  the  basis  of  our  method,  it  is  these 
equations  that  we  will  solve  numerically  in  Section  3.  Note  that  we  can  use 
(2.4)  to  eliminate  0  from  (2.5);  this  leads  to  the  initial-value  problem 
(P* )  for  X. 

Our  procedure  for  solving  (P)  is  based  on  solving  (P*)  for  t  >  0  and 

p  in  the  fixed  interval  [-a, a],  the  support  of  uQ.  (We  will  give  a 

uniqueness  theorem  to  show  that  boundary  conditions  at  p  ■  ta  are  not 

needed.)  Let  X  be  a  solution  of  (P*)>  Since  X  ■  1  at  t  ■  0,  there 

P 


exists  a  T  >  0  such  that 


X  >  0  on  [0,T)  x  [-a, a] 
P 


C±(t)  -  X(t,±a)  , 

-  { ( tf x ) s  C  (t>  <  x  <  C.< t>,  0  <  t  <  T}  . 

Then  for  each  t  e  [0,T)  the  napping  p  X(t,p)  is  a  bijection  of  [~a,a] 
onto  the  interval  tC_ < t) ,C+(t)] ;  letting  P(t,x)  be  such  that 

X(t,P(t,x))  “  x 

(i.e.,  P(t,')  is  the  inverse  of  X(t,»)),  we  define 


u(t,x) 


0(t,P(t,x>>,  (t,x)  e  n 


(2.6) 


^0,  otherwise 

with  0  given  by  (2.4).  He  then  expect  the  resulting  function  u  to  be  the 
weak  solution  of  ( P) .  He  now  show  that  this  expectation  is  indeed 
justified.  To  avoid  repeated  hypotheses  we  assume,  for  the  remainder  of  the 
section,  that 

u0  >  0  on  ("8,8),  u0  ■  0  otherwise,  uQ  e  C(R)  n  C2[-a,a)  .  (2.7) 

Further,  we  will  use  the  term  solution  of  (P  )  on  [0,T]  for  a  solution  X 
on  10, T)  x  [-a, a)  with  X,  Xfc,  Xp,  Xpp,  Xfcp,  and  Xppp  continuous  on 
tO,T)  x  [-a, a] »  if,  in  addition,  Xp  >  0  on  [0,T)  x  (-a, a] ,  then  X  is 

£SSHiiK> 

Theorem  1  (Consistency),  bet  X  be  a  regular  solution  of  (P* )  on  10,T), 
and  let  u  be  defined  by  (2.6).  Them 


(i) 

u  >  0  on  flT,  u  *  0 

otherwise; 

(ii) 

u(x,0)  -  Uq(x)  for 

x  S  Ri 

(iii) 

ut  ■  <«2)xx  ia.  V 

(iv) 

t±(t)  -  -2ux(t,?±(t)) 

and  ±C±(t)  >  0  |o£  t  e  [0,T). 

Proof.  Assertion  (i)  follows  fro*  (2.4),  (2.6),  (2.7),  and  the 
inequality  Xp  >  0>  (ii)  is  a  consequence  of  (2.6),  (2.7),  and  the 
identities  P(0,x)  “  x,  Xp(0, p)  ”  1. 

Next,  <  F  ) ^  and  (2.4)  imply  (2.5)  and  (2.3),  while  (2.1)  is  a  consequence 
of  (2.6).  By  (2.1)  and  (2.5), 


*t  -  -2<ux)  .  Xtp  -  -2(uxx)  Xp  .  (2.8) 

where  we  have  used  the  notation  (2.2).  Also, 

°t  *  (ut>#  +  <ux>\  "  <ut  -  2ttx)#  ' 
and  hence  (2.3)  and  (2.8 )2  imply 

0  *  X  (u  -  2u2  -  2uu  )  ■  X  [u.  -  (u2)  ]  i 

p  t  x  xx  p  t  xx 

since  Xp  >  0,  we  have  (iii). 

Equation  (2.8 )1  leads  to  the  first  of  (iv).  Finally,  since 
u0(±a)  -0,  (P*) i  implies 

C. (t)  -  X  <t,±a)  -  ~2u’ (±a)X  (t,±a)"2  »  (2.9) 

it  up 

but  by  (2.7),  ±u^(±a)  <  0,  and  the  remainder  of  (iv)  follows.  This 

completes  the  proof. 

Theorem  2  (Uniqueness) .  For  any  T  >  0  there  exists  at  most  one  regular 
solution  of  (P  )  on  [0,T). 

Proof,  bet  X  and  Y  be  regular  solutions  on  (0,T).  Choose 
tQ  e  ( 0,T) i  it  suffices  to  show  that  X  -  Y  on  [0,tQ]  x  !“•*«].  By  (P  )^, 

Xt  "  “<uoxp2)p  _  uoxp2'  Yt  ”  "<uoYp2)p  "  uoYp2  * 

Xf  we  subtract  these  equations,  multiply  by  X  -  Y,  and  integrate  fro* 

p  *  -a  to  p  ■  a,  we  obtain,  after  an  integration  by  parts  using 

uQ(±a)  ■  0, 


?  fc  1-.  «-»>'* 


/!.  v*;2  -  <’>«.  -  -  r:2>  ^  • 


(2.10) 
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•iaoe 


<*’2  -  t"2 ) (X  -  Y)  -  -(X  +  T  ><X  -  Y  >W2  <  0 

P  P  P  P  PPPPPP 


tha  f*ret  ttn  on  the  right  a Ida  of  (2.10)  ia  <  0.  Next,  lotting 

♦  -  uXX  +  X  )x‘V2  , 

Op  P  P  P 

wa  find  that 


'/!.  «0<  '  Yp2)(X-y)dP  "  I  /!at(X-Y)23p*dp 

-  -  J  /!.  (x-Y)2*pdp  +  jt(x-x)2t]*a  . 

Sinea  ±u*(±a)  <  0,  we  have  ±*<±a,t)  <  0>  hence 

C(x-Y)2W*  <  o  . 

*  A 

Further,  aa  the  aolutiona  X  and  Y  are  both  regular,  X  and  Y  are 

P  P 

bounded  away  fro*  zero  on  tO,tQ]  *  l-a,a]»  hence  there  la  a  constant  X 
such  that 


on  10, tQ]  x  l -a, a) . 
Let 


'♦,1  <  « 


{<t>  -  /*a  B(t,p)  -  T(t.p»2ap 

and  note  that  5(0 )  -  o,  since  X  and  Y  satisfy  the  same  initial 
condition.  Thus  if  we  integrate  (2.10)  from  0  to  t  and  use  the  above 
remarks,  we  arrive  at  the  Gronwall  inequality 

5(t)  <  K  /J  C(T)dT 

for  0  <  t  <  t0,  which  clearly  implies  X  -  Y  on  10, t^]  *  [-a, a]. 

Remark.  We  have  not  been  able  to  establish  (directly)  existence  for 
Problem  (P*).  We  note  that  existence  can  be  inferred  from  existence  for  (P) 
in  conjunction  with  (1.6). 


Thus  far  all  of  our  results  have  been  local/  as  the  initial  condition 


Xp(°'P)  “  1  only  Insures  >  0  for  sufficiently  small  time*  There  are 
situations  in  which  we  can  give  a  global  result.^ 

Theorem  3.  Suppose  that  u£  <  0  on  [-a, a] .  Let  X  be  a  solution  oi 
( )  on  (0,T).  Then  X  is  regular;  in  fact*  X^  >  1  on  [-a,a]  *  (0,T). 

Proof*  We  let  Z  *  X^  and  differentiate  (P  ) 1  with  respect  to  p  to 
obtain 

3  -12 

Z  «  2unZ  -  6unZ  Z  +  6uiZ  -  2u"Z  .  (2.11) 

tOppOpOpO 

At  t  ■  0  we  have  Z  »  1,  Z„  -  -  0  j  hence  Z.  >  0  at  t  *  0  and  Z  >  1 

p  pp  t 

for  small  t.  Suppose  z  were  ever  equal  to  1.  Then  there  would  be  a  first 
time  t^  e  10, t]  at  which  this  occurs  and  a  point  p^  with  z(tj,pj)  ■  1. 


If  p1  e  (-a 

.a),  we  would  have 

zp  -  0, 

Z  > 

PP 

0/  and  Z^  <  0 

at  (t1/p1). 

which  contradicts 

(2.11)/ 

since 

ttj(pt)  <  0.  if  p1  -  -a. 

then  zt  <  0 

and  Z  >  0  at 

P 

(t^p,)} 

thus. 

since  uQ(-a)  “  0,  u^(-a)  >  0 

and  ujj(-a)  <0,  we  again  contradict  (2.11).  A  similar  argument  applies  at 
x  -  a. 

For  initial  data  uQ  which  is  not  concave  one  cannot  expect  to  have 
*p  >  0  for  all  time.  Knerr  [6],  generalizing  results  of  Aronson  [11],  has 
shown  that  if  Ug(b)  ”  0  for  b  ■  a,  say,  then  the  free  boundary  emanating 
from  a  is  vertical  for  an  interval  0  <  t  <  T  <  •.  At  t  ■  T  there  is  a 
loss  of  smoothness  and  C+ (t)  begins  to  increase.  We  now  show  that  this 
phenomenon  is  related  to  the  vanishing  of  Xp. 


1 

Cf.  Graveleau  and  Jamet  (8],  who  show  that  if  Uq  is  concave/  then  the 
solution  u  of  (P)  will  be  concave  for  all  time.  With  u  <  0,  (2.8) 
clearly  implies  X  >  1. 
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Theorem  4  (Breakdown).  Suppose  that  u^tb)  “  0  and  u|j(b)  >  0  for 
b  “  -a  or  b  ■  a,  and  put 

T  -  [6u£(b))_1  .  (2.12) 

Let  X  be  a  solution  of  ( P  )  on  [0,T).  Then 


(i) 

X(t,b)  =  b 

for  0  <  t  <  T 

(ii) 

xp<t,b)  ♦  o 

as  t  ♦  T. 

If.  in  addition,  X  is  regular  on  [0,T)  and  u  la  defined  by  (2.6),  then 
(ill)  uxx(t,X(t,b))  ♦  -  as  t  T. 

Proof.  Let  z(t)  “  Xp(t,b).  Since  z(0)  -  1,  there  is  a  t0  such 
that  z(t)  >0  for  0  <  t  <  tQ  <  T.  Thus,  by  (2.11),  for  0  <  t  <  tQ, 

z2z*  -  -2u£(b)  (2.13) 

and 

z3(t)  -  1  -  6uJ(b)t  . 

Therefore,  by  continuity,  we  may  take  t0  ■  T  and  (ii)  follows.  Moreover, 

since  z  >  0  on  [0,T),  (P*)  yields  (i).  Finally,  by  (2.8)2  and  (2.13), 

u  (t,X(t,b) )  ■  u"(b)z(t)  3  , 

xx  0 

and  (ii)  implies  (ill). 

Under  the  hypotheses  of  Theorem  3  we  have  the  following  growth  estimate 
2 

for  the  L  norm  of  Xp. 

Theorem  5  (Stability).  Suppose  that  ujj  <  0  on  [-a, a] ,  and  put 

c  ■  8[u£(-a)  -  u^(a)].  Let  X  be  a  solution  of  (P*)  on  [0,T).  Then 

/*  X2(t,p)dp  <  2a  +  ct  (2.14) 

-a  p 

for  0  <  t  <  T. 
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(2.15) 


Proof.  Mote  first  that,  since  uQ(ia)  “  0,  (P  implies 
*t(t,±a)  -  -2u^(±a)x"2(t,±a)  . 

Me  multiply  (P*).  by  X  3X  and  integrate  frost  p  -  -a  to  p  -  a.  After  an 

i  p  pp 

integration  by  parts  and  the  use  of  (2.15)  we  obtain 

-  r  §r  /*  x2dp  -  2 [u'x”1] *  -2  Ja  uftX_3X2  dp  +  2  /*  ul(X_1)  dp  . 

2  dt  -a  p  0  p  -a  -a  0  p  pp  r  '-a  0  p  p  * 

If  we  integrate  the  last  term  on  the  right  by  parts,  we  find  that 

Id  r&  2  r&  *32 

T  —  I  X  dp  +  2  J  u.X  X  dp 
2  dt  '-a  p  J  -a  Op  pp  ^ 


P  '  -a  o  p  pp 


<  4Iul(-a)x“1(t,-a)  -  ui(a)x'1(t,a)] 

up  up 


‘  f  • 


since  u||  <  0,  ±uQ(±a)  <  0,  and  X^  >  1.  the  estimate  (2.14)  follows  upon 
integration . 


Observe  that  one  obtains  also  the  estimate 


/;  /!« *  * 2*  ♦  <*  > 

in  addition,  since  Xp  >  ^  (1.3)  and  (2.15)  yield  the  inequality 

U±(t)|  <  2|u'(±a)|  . 

By  (iv)  of  Theorem  1,  ±X(t,±a)  >  0>  hence  for  each  t  there  is  a  p1 

such  that  X( t ,p^ )  -  0.  The  standard  argvsunt 


|X(t,p)|  -  IP  X  (t,e>dC|  <  [2a  J*  X2(t,C)dC) 
Pj  P  •  P 

and  (2.14)  therefore  yield  the  estimate 


|x(t ,p) |  <  /2a (2a  ♦  ct  . 


1 


Thus  the  width  of  the  support  of  u(*,t)  is  at  most  0(t  ). 
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Remark .  Vasques  [12]  has  shown  that  this  width  is  actually  0(t  )  by 
proving  that  u(t,x)  is  asymptotic  to  a  certain  "fundamental  solution" 

A 

u(t,x)  of  the  form  (3.5).  Such  fundamental  solutions  correspond  to  Dirac 
distributions  at  t  ■  0  and  as  such  the  corresponding  coordinate 

A 

transformation  X( t,p)  is  not  defined.  If,  however,  one  lets  t  denote  the 

ft 

A 

time  for  which  the  support  of  u(t  ,*)  has  width  2a,  and  considers 

ft 

A  A  A 

u  (t,x)  -  u(t  +t,x),  then  the  corresponding  coordinate  transformation  X  is 

ft  ft  ft 

well  defined;  in  fact. 


A  t  +t  3 

-  (— — )  p  . 
a 

A 

Guided  by  the  results  of  [12],  we  conjecture  that  X(t,p)  and  X  (t,p)  are 

ft 

asymptotic  as  t  ♦  •. 

Problem  (P* )  can  be  given  a  weak  formulation,  which  we  now  deduce.  We 


t 

begin  by  writing  (P  ) 1  in  the  form 


Xk  -  -<u.X_2>  -  u'X~2 

t  0  p  p  Op 


If  we  multiply  this  equation  and  (P* )2  by  an  arbitrary  C1  function  w(p)  and 
integrate  from  p  -  -a  to  p  ■  +a  we  obtain,  after  an  integration  of  parts. 


/!.<« t  *  “{’'p2**  -  Vp2*’14'’  ’  °  ■ 

(2.16) 

J*a [X(0,p)  -  p] w( p) dp  -  0  . 

Equations  (2.16)  constitute  the  desired  weak  formulation  of  (P*). 

The  weak  form  (2.16)  admits  an  approximate  formulation  in  terms  of  finite 
elements.  In  this  connection  it  is  important  to  note  that  the  absence  of 
boundary  conditions  allows  one  too  operate  in  the  space  H 1 ( -a , a ) .  We  expect 
these  observations  to  be  of  great  value  in  the  extension  to  higher  dimensions. 
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Remarks 


1.  The  p  ”  constant  trajectories  t  **  x(t,p)  are  analogs  of 
characteristic  curves ,  as  the  degeneracy  in  the  porous  media  actuation  (P ) 1 
propagates  along  such  trajectories.  Indeed,  this  equation  degenerates  at 
u  ■  0,  and 


u( t ,x ( t ,p) )  *  0  if  and  only  if  Ug(p)  "0  , 

at  least  when  the  solution  is  regular  (cf.  (2.4)). 


2.  Problem  (P)  is  a  special  case  of  the  more  general  diffusion  problem 

(2.17) 


“t  ■  *l“'u,,x  ' 


u(x,0)  -  uQ(x)  # 


u,C)  | 


with  (2.16)  degenerate  in  the  sense  that 

0  for  u  >  0 
0  for  u  *  0  . 

As  before,  we  seek  curves  x  -  X(t,p)  along  which  this  degeneracy  propagates. 

Formally,  such  curves  are  generated,  via  (1.6),  by  writing  (2.17)^  as  a  mass 

balance  law  (1.4),  since  this  law  has  the  Lagrangian  form 

u( t,X(t ,p) )  -  X  1(t,p)uft(p) 

P  0 

(cf.  (1.7)).  if  we  do  this,  we  find  that  the  resulting  initial-value  problem 
for  X(t,p)  is: 

*  u_  ,  u. 

"  .i-  {£)  )  . 

(2.18) 


0  p  p  p  p 


X(0,p)  -  p  . 

As  an  example,  the  porous  media  equation  is  often  considered  in  the  form 

u  ■  <«")_  (■  >  2)  » 

t  xx 

for  this  equation  (2.18)^  becomes 


m-2 


■u.  u. 

_ _ 2— 

t  _n-1  > 


P 


P  P 
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3.  numerical  solutions 


Here  m  describe  a  simple  difference  scheme  for  the  initial-value  problem 
(P*).  (We  will  actually  wort  with  (2.4),  (2.5)^  rather  than  (P  )j.)  For 
convenience,  we  assume  that  u0  is  symmetrict  uQ(x)  -  u0(-x).  Then  by 
symmetry  we  can  restrict  ourselves  to  0  <  p  <  a  provided  we  impose  the 
additional  conditions 


X(0,t)  5  0,  U(0, t)  s  0 
P 


(3.1) 


We  want  to  solve  the  equations 


-  -x-1u  , 

t  P  P 


0  "  X  *u_ 
P  0 


(3.2) 


We  choose  At  >  0  and  Ap  -  aN  1  for  seem  integer  N,  introduce  the 

mesh  t,  -  iAt,  i  -  0,1,2,...,  p  -  nAp,  n  -  0,1, ...,N,  and  write  f*  for 
in  n 

the  value  of  a  function  f  at  the  mesh  point  (tirpn).  We  approximate  (3.2) 
by  the  difference  scheme: 


x^-x1  t  o1 

n _  n  n 


.i-M  _  n 


MX1'  _MXi+1  * 

n  n 


(3.3) 


In  these  formulae  -  uQ(pn),  while  1^  and  ^  represent  the  spatial 
difference  operators 


Lnf  "  2Ap 


^n+1  *n-1 ' 


fN-2  “  4fH-1 ' 


N  f 
n 


1 

2Ap 


0, 


f n+1  *n-1' 


4Vl  +  fN-2' 


2f1' 


1  <  n  <  N  -  1 

n  ■  H 

n  *  0 

1  <  n  <  N  -  1 
n  ■  W 
n  •  0  . 


0.4) 
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2 

The  formulae  <3.4 )  are  accurate  to  0(Ap  )  for  functions  satisfying 

(3.1)  with  U(t,a)  s  0.  since  the  method  (3.3)  is  partially  explicit,  we 

2  2 
would  expect  it  best  to  take  At  *  (Ap)  »  with  this  choice  we  expect  0(Ap  ) 

accuracy. 

The  following  explicit  solution  to  the  initial-value  problem  (P)  -  for 
initial  data  a  Dirac  distribution  -  is  given  by  Pattle  [13] » 

2 

(  Mt^M  -  r),  lyl  <  X(t> 

;<«.t)  -  {  x<t) 

^  0,  lyl  >  X(t)  (3.5) 

1_ 

y  -  x/xQ,  X(t)  -  <t/t0)3,  xQ  -  r(5/2)//r  -  “ ,  tQ  -  x*/12  . 

Jk 

We  attempted  to  approximate  the  solution  u(t,x)  ”  u(  t+1 , x) .  Thus  we  have 

a  “  (12xq)1^3  ■  91//3.  The  results  are  presented  in  Tables  1  and  2. 

Xn  Table  1  we  give  the  approximate  and  theoretical  values  of  C.(t), 

▼ 

together  with  the  relative  errors,  for  a  sequence  of  times.  These 

calculations  are  performed  with  ten  subdivisions,  so  that  Ap  ■  0.208.  we 

observe  that  even  with  this  fairly  crude  mesh  the  free  boundary  is  tracked 

quite  accurately.  In  Table  2  we  present  one  of  our  calculations  to  determine 

2 

the  rate  of  convergence.  The  results  confirm  the  expected  rate  of  (Ap)  . 


-14- 


Table  1 


A  T 

£  (t)  end  £  (t),  the  approximate  and  theoretical  positions  of  the  right- 

2 

hand  free  boundary  for  dp  ■  0.208,  At  “  <Ap)  . 


t 

CA(t) 

T 

£  (t> 

relative  error 

5.327 

3.63826 

3.63274 

0.00152 

9.653 

4.43406 

4.42904 

0.00113 

13.980 

5.01550 

5.01093 

0.00091 

18.307 

5.48638 

5.48218 

0.00079 

22.634 

5.88782 

5.88391 

0.00063 

26.960 

6.24088 

6.23721 

0.00059 

31.287 

6.55793 

6.55445 

0.00053 

35.614 

6.84696 

6.84365 

0.00048 

39.941 

7.11343 

7.11027 

0.00044 

44.267 

7.36131 

7.35827 

0.00041 

Table  2 

A  T 

C  (1.108)  and  £  (1.108),  the  approximate  and  theoretical  positions  of  the 

2 

right-hand  free  boundary  at  t  ■  1.108  for  At  “  (Ap)  . 


Ap 


£*< 1.108) 


£T ( 1 . 08 ) 


error 


0.1040042 

0.0520021 

0.0260011 

0.0130005 


2.152773 

2.152591 

2.152546 

2.152535 


2.152531 

2.152531 

2.152531 

2.152531 


0.000242 

0.000060 

0.000015 

0.000004 


Problem  (P)  with 

f  2  wx  .  t 
j  cos  r-,  x  <  1 

u0<x>  (3.6) 

V  0,  otherwise 

is  discussed  by  Aronson  t 1 1 1 »  who  shows  that  the  corresponding  free  boundaries 
are  vertical  for  an  interval  0  <  t  <  T,  a  result  consistent  with  the 
conclusions  of  Theorem  4,  where  for  (3.6), 


15 


T  -  (3*2)"1  *  .03377 


We  also  performed  a  numerical  experiment  for  the  initial  data  (3.6)f  the 
results  are  shown  in  Table  3.  It  is  seen  that  the  free  boundary  is  indeed 
roughly  vertical  for  t  <  T. 


Table  3 

CA(t),  the  approximate  position  of  the  right-hand  free  boundary  for 
Uq(x)  -  cos2  Here  Ap  -  0.025,  At  ■  (Ap)2. 

t  cV> 

0.00375  1.0000016 

0.99750  1.0000039 

0.01500  1.0000130 

0.02000  1.0000293 

0.02500  1.0000752 

0.03000  1.0002330 

0.03125  1.0003195 

0.03250  1.0004445 

0.03375  1.0006275 

0.03500  1.0008996 

0.03625  1.0013092 

0.03750  1.0019303 


] 

I 


4.  lx tension  to  tf*. 

The  porous  media  problem  in 
u(t,g),  t  >  0,  j  8  l",  such  that 


*”  consists  in  finding  a  scalar  function 


-  A  (u2) 
2 


«V 


«<°*2>  “  V**  ' 

(Bars  and  in  what  follows  A,  div,  and  7,  respectively,  dsnots  the 
Laplacian,  divergence,  and  gradient  operators  in  h". ) 

Guided  by  our  one-dimensional  analysis,  we  rewrite  (P_) -  as  a  balance 

ll  1 


law 


with 


u.  ♦  div  (uv)  »  0 
t  x  *■ 


v  “  -27  u  , 

~  x  ' 

aw 

and  we  take  x(t,j)  to  be  the  solution  of  the  initial-value  problem 

Vt'B>  “  *<°'£>  •  £  • 

Procedi ng  as  before,  we  define 

0(t,p)  «  u(t,x(t,g)) 


and  note  that 


t<et 


(ut  + 


7  u) 
x 


-(u  div  v) 

xt* 


(4.1) 


(4.2) 


and  assume  that  det  Z  >  0.  Then,  using  the  identities  (ef.,  e.g.,  [10J,  p. 
77) 


(det  Jg)t  -  (det  Z)tr(gtg’1)  , 
tr^-1)  -  <divx  v)* 

in  conjunction  with  (4.2),  we  conclude  that 

(O  det  pt  -  0 
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and  hence  that 


U 


(det  Z)  \ir 

L 


On  the  other  hand,  by  the  chain-rule, 

V  O  -  ZT(V  u)*  . 

B  ~  5 

Thus,  using  (4.1),  we  arrive  at  the  following  initial-value  problem  for  X: 

*  -M'  «  "  VX)_1u  , 

<Pn> 

X(0,g)  -  g  > 


where  V  «  V  . 

E 

* 

A  careful  analysis  of  (P  )  is  beyond  the  scope  of  this  paper.  Our 

n 

ultimate  hope  is  to  show  that  when 

Uq  >  0  on  A,  u0  -  0  otherwise,  (4.3) 

with  A  compact  and  connected.  Problem  (PR)  reduces  to  solving  (P^)  on  the 
fixed  region  A  for  all  time. 

* 

In  view  of  (4.3),  ufl(g)  •  0  for  Jg  €  3A.  Thus,  by  lPn),  £t(t,g)  =  0 
at  any  g  e  3a  for  which  VuQ(g)  **  0,  at  least  as  long  as  the  solution 
remains  regular. 

Remarks. 

1.  The  p  “  constant  curves  t  **■  X(t,g)  are  analogs  of  bicharacteristic 
curves  in  the  theory  of  partial  differential  equations  (cf.  Remark  1  at  the 
end  of  Section  2). 

2.  For  Problem  (Pn)  with  n  >  2,  regularity  of  the  free  boundary  f  is 
essentially  an  open  question  (cf.  [14]).  The  formulation  (P^)  might  be 
useful  in  attacking  this  question,  as 

F  -  (X(t,g)  i  g  «  JA,  t  >  0}  . 

Moreover,  since  regularity  is  a  local  question,  the  problem  of  proving  that 
X(t,*)  is  a  bijection  is  trivial!  it  follows,  at  least  locally,  from 
X(0,g)  -  g,  where  we  have  chosen  the  time  scale  with  t  ■  0  the  time  near 
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which  regularity  is  sought. 

3.  For  the  sore  general  equation 

u  ■  A  (u*>  {«  >  2) 

• 

the  first  of  (P  )  is  replaced  by 
n 

vxTx  -  -■u**“2vu  . 
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